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The spin- 1/2 quantum antiferromagnet on the Kagome lattice provides a quintessential example 
in the strongly correlated electron physics where both effects of geometric frustration and quantum 
fluctuation are pushed to their limit. Among possible non-magnetic ground states, the valence bond 
solid (VBS) with a 36-site unit cell is one of the most promising candidates. A natural theoretical 
framework for the analysis of such VBS order is to consider quantum states on a bond connecting 
the nearest-neighboring sites as fundamental quantum modes of the system and treat them as 
effectively independent "bond particles." While correctly describing the VBS order in the ground 
state, this approach, known as the bond operator theory, significantly overestimates the lowest spin 
excitation energy. To overcome this problem, we take a next logical step in this paper to improve 
the bond operator theory and consider extended spin clusters as fundamental building blocks of 
the system. Depending on two possible configurations of the VBS order, various spin clusters are 
considered: (i) in the VBS order with staggered hexagonal resonance, we consider one spin cluster 
for a David star and two spin clusters with each composed of a perfect hexagon and three attached 
dimers, and (ii) in the VBS order with uniform hexagonal resonance, one spin cluster composed of a 
David star and three attached dimers. It is shown that the majority of low-energy spin excitations 
are nearly or perfectly fiat in energy. With most of its weight coming from the David star, the 
lowest spin excitation has a gap much lower than the previous value obtained by the bond operator 
theory, narrowing the difference against exact diagonalization results. Despite the better agreement, 
however, the fact that the lowest spin excitation mostly comes from the David star structure is at 
odds with the bond operator theory as well as previous series expansion studies. Expecting that the 
true low-energy spin excitations have contributions from both the David star and its surroundings, 
we propose to use two super spin clusters in future work treating spin excitations from the both 
contributions on a more equal footing. 



I. INTRODUCTION 

Quantum antiferromagnet in two spatial dimensions 
has been a subject of intense research for multiple rea- 
sons. On one hand, the antiferromagnetic order param- 
eter itself does not commute with the antiferromagnetic 
exchange interaction in the Hamiltonian so that there ex- 
ist inherent quantum fluctuations in the system. On the 
other hand, effects of quantum fluctuations are amplified 
as the size of the spin moment becomes smaller and the 
spatial dimension of the system gets lower. The antifer- 
romagnetic spin- 1/2 chain is the limiting situation where 
both reach their respective minimum possible values. In- 
deed, the antiferromagnetic spin- 1/2 chain demonstrates 
the strongest effect of the quantum fluctuation, which 
makes the ground state a gapless spin-liquid/^!^ While 
interesting in many aspects, the one-dimensional quan- 
tum antiferromagnet is too much susceptible to quantum 
fluctuations that it hardly displays any long-range order. 

In this context, the two-dimensional quantum anti- 
ferromagnet gives rise to perhaps the most interesting 
regime where the strength of quantum fluctuation is just 
so right that there is a delicate balance between order 
and disorder. On simple lattices, magnetically ordered 
states are often found to be energetically favorable. For 
example, the two-dimensional quantum antiferromagnet 
becomes ordered in the square lattice with staggered spin 



moments. The ground state in this case is known as the 
Neel state. Realizing quantum disordered states in two 
spatial dimensions requires a delicate way of tilting the 
balance between order and disorder. Geometric frustra- 
tion is often envisioned as an effective way to do just 
that. 

The Kagome lattice is regarded as one of the most 
geometrically frustrated systems in two dimensions ow- 
ing to the corner-sharing triangular structure, and there- 
fore expected to be the best place to look for the quan- 
tum disordered state in spatial dimensions higher than 
one. Indeed, in a recent experimental discovery of 
an "ideal" Kagome lattice realized in Herbertsmithites, 
ZnCu3(OII)6Cl2,— it is shown that the material remains 
magnetically disordered down to 50 mK while its Curie- 
Weiss temperature, a measure of the spin exchange en- 
ergy, is about 300 K. This experimental finding strongly 
suggests that the ground state might be spin disordered 
even at zero temperature, providing a concrete example 
of fully quantum disordered states in two spatial dimen- 
sions. Despite this strong suggestion for the existence, 
however, the precise nature of the spin disordered state 
is still unknown. 

To construct a precise theory for the spin disordered 
ground state in the Kagome lattice antiferromagnet, it 
is crucial to use lessons obtained from numerical studies. 
Following are some of the most important lessons from 
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FIG. 1. (Color online) Two spin-singlet covering configurations of the valence bond solid (VBS) ground state with a 36-site 
unit cell. The left figure shows the VBS with staggered hexagonal resonance and the right figure shows the one with uniform 
hexagonal resonance. Each dimer, denoted as a thick line in the figure, indicates the position of a spin-singlet bond. The unit 
cell with 36 sites (or 18 dimers) is indicated by the enclosing dashed line. Note that the two VBS order states have different 
symmetries. With respect to the rotation around the center of the pinwheel structure, the VBS order with staggered hexagonal 
resonance has a sixfold rotation symmetry while that with uniform hexagonal resonance is threefold. 



exact diagonalization<^"— (i) The ground state is mag- 
netically disordered, (ii) The spin-0 excitation seems to 
have a very small gap or is completely gapless. (iii) The 
spin-1 excitation has a small, but finite gap, which is 
0.164 J in the = 36 system and estimated to be around 
0.05 — O.IJ in the thermodynamic limit. Here, J is the 
spin exchange energy. 

There are two possible candidates for the spin dis- 
ordered ground state. ^"^"^ One is the spin liquid state 
where spin moments are disordered without breaking any 
lattice translation symmetryJ^""— There are two varia- 
tions within the spin liquid scenario, difi^ering in whether 
the spin excitation is gapped or gapless. Considering 
the results from exact diagonalization and recent den- 
sity matrix renormalization group (DMRG) studiea^^*^, 
a gapped version of the spin liquid has a potential to be 
the true ground state. The other candidate is the so- 
called valence bond solid (VBS) order where spins are 
paired up with one of its nearest neighbors to form a 
spin singlet, consequently breaking the lattice transla- 
tion symmetryJ^"— While earlier exact diagonalization 
as well as series expansion studies favor the VBS order, 
a recent study using the DMRG method^^ seems to indi- 
cate that the spin liquid may be closer to the true ground 
state. 

Resolving the issue as to which state is really the true 
ground state for the pure Heisenberg model is highly 
non-trivial and numerically subtle since many competing 
ground states are very close in energy. Setting aside this 
issue, however, what is important is the fact that, due 
to the very proximity between various competing ground 
states, even a small deviation from the pure Heisenberg 
model Hamiltonian can cause significant changes in the 
ground state property. For example, it is shown in the 
same DMRG study mentioned above^^ that a slight re- 



duction in J by less than 2% in a particular pattern can 
in fact stabilize the VBS order with a 36-site unit cell 
rather than the spin liquid. In addition, there are various 
important experimental factors that modify the Hamilto- 
nian and thus play crucial roles in determining the true 
ground state. Such factors include the Dzyaloshinskii- 
Moriya interaction, the lattice deformation, the presence 
of impurity spins, and so on. In this context, it is con- 
structive to conduct a comprehensive investigation of all 
possible spin disordered states and study their properties. 

In this paper, we focus on the VBS state with a 36- 
site unit cell illustrated in Fig. [T] Spin singlet bonds, 
denoted as thick lines in the figure, represent the VBS 
ordering pattern with the 36-site unit cell; this is be- 
lieved to be the most favorable configuration by many 
of the previous studies ji^^—ii^ Schematically speaking, 
the most favorable VBS order is the one with the hon- 
eycomb structure of perfect hexagons which in turn sur- 
round the pinwheel structure in the David-star shaped 
region at the core. Here the perfect hexagon means that 
the sides of a hexagon is alternatively covered with three 
spin-singlet pairs. Fig. [T] shows two possible honeycomb 
arrangements of the perfect hexagons with one having the 
staggered hexagonal resonance and the other the uniform 
hexagonal resonance. In the VBS order with staggered 
hexagonal resonance, perfect hexagons have an alternat- 
ing dimer covering pattern (denoted as filled and empty 
circles inside the hexagon) encircling the pinwheel struc- 
ture. Meanwhile, in the VBS order with uniform hexago- 
nal resonance, all perfect hexagons have the same dimer 
covering pattern. Note that the pinwheel structure it- 
self can have two different orientations, which have ex- 
actly the same energy. The degeneracy between two such 
states is not lifted even after the pinwheel is coupled with 
surrounding dimers. This is in fact the origin of the gap- 
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less spin-0 excitations in the VBS scenario. We consider 
both the VBS order with staggered and uniform hexag- 
onal resonance in this paper since they are equally good 
candidates for the most stable VBS order. 

The VBS state with uniform hexagonal resonance was 
previously studied by some of the authors within the 
bond operator theory (BOT).^^ The bond operator the- 
ory provides a natural theoretical framework for the anal- 
ysis of the valence bond solid state by regarding quantum 
states on a bond connecting the nearest-neighboring sites 
as fundamental quantum modes of the system. At the 
mean-field level, such quantum modes are treated as ef- 
fectively independent "bond particles." While correctly 
producing the VBS order in the ground state, however, 
the bond operator theory turns out to be quantitatively 
inaccurate. In particular, the spin-1 excitation gap is 
predicted to be much larger than that obtained from ex- 
act diagonalization of the 36-site cluster while the ground 
state energy is not too bad considering the level of ap- 
proximation. 

In addition, a peculiar property of the bond operator 
theory is that the lowest spin-1 excitations entirely origi- 
nate from the surroundings [light blue (gray) in the right 
figure of Fig. [1] that encircle the core pinwheel struc- 
ture. This is a consequence of the bond operator theory 
that core dimers [dark blue (gray) in the right figure of 
Fig. [T] are completely decoupled from the surrounding 
ones. While this conclusion may be fully consistent with 
series expansion results, the large discrepancy in the low- 
est spin-excitation gap energy leaves the accuracy of the 
bond operator analysis in question. Moreover, the bond 
operator theory predicts that the spin-1 excitations from 
the pinwheel structure are all degenerate with the energy 
gap being exactly IJ. This is clearly not correct. 

While the problem certainly arises from the coarse na- 
ture of the mean-field approximation, there is no system- 
atic way to improve the mean-field theory, especially for 
the "hard-core" constraint imposing that no two bond 
particles can occupy the same dimer site. In this pa- 
per, instead of looking for a perturbative way to include 
higher-order diagrams, we set out to improve the bond 
operator theory by extending the fundamental building 
block of the system from a bond to a more extended clus- 
ter of spins. 

To determine which spin clusters are most appropriate, 
it is worthwhile to investigate the nature of the inter- 
dimer interaction carefully. To this end, it is important 
to note that all inter-dimer interactions can be classi- 
fied into two basic units according to their geometric 
structures: (i) topologically orthogonal and (ii) empty 
triangle structure. See Fig. [2] for illustration. Extended 
spin clusters can be constructed by the combination of 
these two basic units. Next, by using exact diagonal- 
ization, we investigate how each basic structure affects 
the low-energy spectrum of the extended spin clusters. 
As a result, it is shown that the empty triangle struc- 
ture mainly contributes to the lowering of the ground 
state energy while the topologically orthogonal structure 



stabilizes the ground state itself by reducing quantum 
fluctuations about the VBS order. 

We first test the bond operator theory as to how well it 
captures the effects of the topologically orthogonal as well 
as the empty triangle structure. The test shows that the 
bond operator mean-field theory is quite accurate for cap- 
turing the effects of the empty triangle structure while it 
is not so in the case of the topologically orthogonal struc- 
ture. Based on this result, we conclude that the most 
appropriate spin cluster is the one containing the max- 
imal number of the topologically orthogonal structures 
allowed by the VBS order. An important assumption 
here is that, similar to the bond operator mean-field the- 
ory, the effects of all remaining empty triangle structures 
are well captured by the mean-field theory formulated in 
terms of the quantum modes formed in the extended spin 
clusters. For convenience, let us call this approach the 
spin cluster operator theory (SCOT). 

Figure [3] shows spin cluster coverings for each type of 
the VBS order depicted in Fig. [1] As one can see, the 
VBS order with staggered hexagonal resonance is covered 
with one spin cluster for a David-star structure (denoted 
as D) and two spin clusters with each composed of a per- 
fect hexagon and three attached dimers (denoted as HI 
and 112). Meanwhile, the VBS order with uniform hexag- 
onal resonance is covered with one spin cluster composed 
of a David-star and three attached dimers (denoted as C) 
and nine dimers for spin-singlets (denoted as SI, • • • , S9). 
Here spin clusters are chosen according to the following 
three conditions: (i) the VBS order formed inside each 
individual spin cluster should be consistent with that of 
the whole lattice, (ii) the shape of spin clusters should 
preserve the symmetry of the VBS order, and (iii) spin 
clusters should include as many inter-dimer interactions 
with the topologically orthogonal structure as possible 
within themselves. There is one additional practical con- 
straint in the choice of the above spin clusters, which 
comes from the technical problem that the size of the 
Hilbert space grows exponentially as the number of spins 
increases. For this reason, in this study, we choose an 
appropriate subset of the unit cell as a spin cluster. For 
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FIG. 2. (Color online) Two diflerent structures of the inter- 
dimer interaction: (i) topologically orthogonal and (ii) empty 
triangle structure. Alphabetical and numerical indices denote 
the individual spin and the dimer, respectively. 
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FIG. 3. Spin cluster covering for the VBS order with staggered hexagonal resonance (left figure) and that with uniform 
hexagonal resonance (right figure). The unit cell of the VBS order with staggered hexagonal resonance consists of three spin 
clusters that include a David star, D, and two extended hexagons, HI and H2. Similarly, the unit cell of the VBS order with 
uniform hexagonal resonance is composed of a core cluster, C, and nine surrounding dimers, SI - S9. 



future work, we propose to investigate the "super clus- 
ter" that contains all topologically orthogonal structures 
within the whole unit cell. Details of our proposal are 
discussed later in this paper. 

The ground state energy as well as the spin-1 exci- 
tation energy spectrum can be computed via mean-field 
analysis of the spin cluster operator theory. For the VBS 
order with staggered hexagonal resonance, the ground 
state energy is given by — 0.428J per spin which is not 
too bad in comparison with the 36-site exact diagonal- 
ization result, —0.438 J, considering the coarse nature of 
the mean- field theory. The spin-1 excitation spectrum is 
found to be almost flat with very narrow band widths 
in the low energy sector. The lowest spin-1 excitation 
occurs at the K point in the first Brillouin zone with the 
energy gap being 0.323J. Compared to the mean-field 
result of the bond operator theory for the same VBS or- 
der, 0.833J, the above spin-1 excitation energy gap is in 
much better agreement with exact diagonalization. For 
the VBS order with uniform hexagonal resonance, the 
core cluster, denoted as C in Fig. [3l is completely decou- 
pled from the surrounding dimers, resulting in exactly 
the same ground state energy as that of the bond oper- 
ator theory, — 0.422J per site. Spin-1 excitations consist 
of two entirely separate groups with one coming from the 
core cluster and the other from the surrounding dimers. 
The lowest spin-1 excitation occurs from the core clus- 
ter with its energy gap being 0.243 J. This energy gap 
is also much reduced from the BOT result for the same 
VBS order, 0.776J. 



Despite these improvements, however, the spin-1 exci- 
tation energy gap is still somewhat higher than the exact 
diagonalization result. Also, an important feature of the 
low energy excitations is that they predominantly origi- 
nate from the triplet states within the pinwheel structure 
in Fig. 131 This fact is at odds with the bond operator the- 
ory as well as previous series expansion studies. In com- 
bination with the reduced, but still larger size of the spin 
gap, this suggests that the effects of the topologically or- 
thogonal inter-dimer interaction structures are not fully 
incorporated. Expecting that the true low-energy spin 
excitations have contributions from both the pinwheel 
and its surroundings, wc propose to use super spin clus- 
ters in future work, which treat spin excitations from 
the both contributions on a more equal footing. The 
proposed super clusters are chosen so that they include 
all topologically orthogonal inter-dimer structures within 
the 36-site unit cell. Details of the proposal are discussed 
in the Discussion section. 

The rest of the paper is organized as follows. In 
Sec. ini we study physical effects of the inter-dimer 
interaction in two basic units with different geometric 
structures: (i) topologically orthogonal and (ii) empty 
triangle structures. The bond operator theory is tested 
on such basic unit structures in order to clarify how well 
it captures the effects of each inter-dimer interaction. 
The spin cluster operator theory is proposed as a 
solution to the problems emerging in the bond operator 
theory. In Sec. IIIII and llVi the formalism for the spin 
cluster operator theory is constructed for both the VBS 
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order with staggered and uniform hexagonal resonance. 
Resuhs of the spin cluster operator theory are found 
in Sec. |V] and I VII In Sec. IVIIi we propose two super 
clusters of spins; the analysis of which can further 
improve the accuracy of the spin cluster operator theory. 
Finally, we conclude in Sec. IVIIII 



II. MOTIVATION 
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As mentioned in the Introduction, in this paper, the va- 
lence bond solid state with a 36-site unit cell is regarded 
as a candidate ground state for the Kagome-lattice anti- 
ferromagnetic Heisenberg model (KLAHM): 

i/^jJ^S.-S,, (1) 

where denotes that i and j are the nearest neigh- 

bors. Here, the spin exchange energy, J, is positive since 
we are interested in the antiferromagnetic model. The 
VBS order with a 36-site unit cell has a dimer cover- 
ing pattern represented by a honeycomb array of perfect 
hexagons with a pinwheel structure at the center of each 
cell. This particular pattern of the VBS order is generally 
believed to be the most favorable among various possi- 
bilities while the underlying physical reason why this is 
so varies from theory to theory. 

With the VBS order composed of the two structures, 
(i) a honeycomb array of perfect hexagons and (ii) a pin- 
wheel structure at the center of each honeycomb cell, dif- 
ferent viewpoints occur regarding which structure plays 
a more important role for the stability of the VBS order. 
By using a duality mapping to the fully frustrated Ising 
model, Nikolic and Senthil^^ advocated a viewpoint that 
the most preferred dimer covering is the one simultane- 
ously maximizing the number of perfect hexagons and 
pinw heels (or the 12-bond star-shaped flippable loops as 
called in the paper). Another viewpoint is provided by 
Singh and Huseii who, by using series expansion, argue 
that the ground state energy is minimized when every 
three empty triangles are bound into perfect hexagons. 
Such formed perfect hexagons are furthermore arranged 
in such a way that empty triangles between the nearest 
perfect hexagons share as many dimer bonds as possible. 
The resulting dimer covering is exactly the same honey- 
comb array of perfect hexagons as before. In this picture, 
the pinwheel structure at the center of the unit cell is 
just a simple byproduct of the overall covering pattern. 
Some of the authors in this paper proposed a yet differ- 
ent viewpoint that the number of topologically orthogo- 
nal inter-dimer structures is maximized for the stability 
of the dimer covering pattern. Both the honeycomb ar- 
ray of perfect hexagons and the pinwheel structure at the 



TABLE I. (Color online) Two basic units for the inter-dimer 
interaction: (Tl) topologically orthogonal and (El) empty 
triangle structure. Thick (blue) line segments denote where 
spin singlet dimers are located. The table shows exact diago- 
nalizaion results for the ground state energy per spin, tgr/ J, 
the lowest spin-1 excitation energy gap, As^i/J, and the 
number of states within IJ from the ground state in energy, 
Nij. Numbers near thick and thin lines in the figure denote 
the spin-spin correlation value, (Si ■ Sj), for the ground state. 
Corresponding numbers for the lowest spin-1 excitation are 
given in parenthesis. 



center of the unit cell are the direct consequence of such 
maximization. In this paper, we take this approach. 

To reveal why this approach is physically appealing, 
let us investigate how the two different types of the basic 
inter-dimer interaction unit, (i) topologically orthogonal 
and (ii) empty triangle structure, affect the ground state 
as well as the low-energy excitations. We first show ex- 
act diagonalization results for each basic inter-dimer in- 
teraction unit in Table HI Then, the results for further 
complicated spin clusters extending the two basic units 
are presented in Table [Tl] and IIIII 

The most important feature of the spin-singlet ground 
state in the topologically orthogonal unit structure, de- 
noted as Tl, is that the exact ground state is a simple 
direct product of the two spin-singlet states formed in 
each dimer bond depicted as thick (blue) lines in Table ID 
Considering that such direct product state is simply the 
ground state with the inter-dimer interaction completely 
turned off, it means that there is absolutely no quantum 
fluctuation about the valence bond solid order. In fact, 
no matter how large the spin cluster may become, the di- 
rect product state of spin singlets can be shown to be an 
exact energy eigenstates of the Heisenberg-model Hamil- 
tonian (while not necessarily the ground state) so long 
as all inter-dimer connections are topologically orthogo- 
nal. The Shastry-Sutherland lattice provides a realistic 
example where actually all inter-dimer connections are 
topologically orthogonal in the thermodynamic limitj^ 
In this situation, it is known that the VBS state re- 
mains as the exact ground state so long as the inter-dimer 
spin exchange energy is less than a certain critical per- 
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centage (estimated to be around 0.7) of the intra-dimer 
counterpart It is worthwhile to mention that an ac- 
tual material, SrCu2 (603)2, has been discovered to be a 
faithful realization of the Shastry-Sutherland lattice with 
the VBS state being the ground state. Coming back to 
the Kagome lattice problem, what is important here is 
that the VBS order can be stable in the topologically or- 
thogonal inter-dimer interaction structures owing to the 
absence of the quantum fluctuation. 

In the case of the empty triangle unit structure, de- 
noted as El in Table HI the ground state wave function is 
no longer a direct product of spin singlet wave functions. 
The overall characteristics of the ground state is, how- 
ever, still quite well captured by the VBS order with three 
spin singlets formed in each dangling bond attached to 
the empty triangle. Here, we use a simple criterion that 
the closer the spin-spin correlation value, (Sj • Sj), is to 
-0.75 (0.25), the better the state is described by the spin 
singlet (triplet) wave function. Meanwhile, the ground 
state energy per spin in the empty triangle unit structure, 
— 0.417J, is much lower than that of the topologically or- 
thogonal counterpart, —0.375 J. In combination with the 
loss of the exact VBS order, this fact can be interpreted 
as if the empty triangle structure allows more quantum 
fluctuations in exchange for a lowered ground state en- 
ergy. So far, the conclusion is that the topologically or- 
thogonal structure stabilizes the VBS order by removing 
quantum fluctuations while the empty triangle structure 
tries to lower the ground energy by allowing some. It is 
shown below that this conclusion remains valid for spin 
clusters larger than the two basic structures studied here. 

Table |TT] shows a variety of spin clusters composed of 
the topologically orthogonal, T2-T6, and the empty tri- 
angle, E2-E5, inter-dimer interaction structures. In spin 
clusters T2-T6 where entire dimer connections are topo- 
logically orthogonal, the ground state wave function is 
exactly given by the direct product of spin-singlet wave 
functions formed in each dimer. Similar to the topolog- 
ically orthogonal unit structure, Tl, there is absolutely 
no quantum fluctuation in the ground state for the VBS 
order so long as all inter-dimer connections are topologi- 
cally orthogonal. This plays a very important role for the 
stability of the VBS order in the Kagome-lattice antifer- 
romagnet. We believe that, given the huge degeneracy 
in the spin configuration of the (semi-classical) ground 
state in the KLAHM, the ground state tries to lift as 
much degeneracy as possible by minimizing the quan- 
tum fluctuation. One way to achieve this goal is to break 
the lattice translational symmetry with a 36-site unit cell 
where the number of topologically orthogonal inter-dimer 
structures is maximized. In view of an effective statis- 
tical model obtained via quantum-to-classical mapping, 
we envision that the effective temperature is less than 
the critical value so that the system tries to minimize 
quantum fluctuations. This viewpoint is similar to that 
of Nikohc and Senthil 

Despite the formation of a pure VBS order, the quan- 
tum fluctuation is not actually completely gone due to 



a two-fold degeneracy in the pinwheel structure. Note 
that all topologically orthogonal spin clusters, (T2)-(T6), 
have a pinwheel-shaped dimer covering structure in the 
middle. Each VBS ground state shown in the Table [ll] 
has a degenerate energy partner that can be obtained by 
reversing the orientation of the dimer covering pattern in 
the pinwheel structure. This fact is actually very impor- 
tant in view of the numerical result that the spin-0 exci- 
tation is known to be either completely gapless or gapped 
with a very small gap. Reversing the dimer covering ori- 
entation in the pinwheel structure changes neither the 
spin value nor the energy, which, therefore, becomes the 
most natural candidate for the spin-0 excitation mode. 

Now, let us switch gears to spin excitations. As one 
can see from the lower left flgure in Table [TTl the spin-1 
excitation gap decreases as the number of topologically 
orthogonal inter-dimer interactions increases from Tl to 
T6. At the same time, more and more spin-1 excitations 
are available within IJ in energy from the ground state, 
leading to a dense low-energy spectrum. It is shown later 
that it is this decrease (as well as the dense low-energy 
spectrum) that the bond operator theory misses, eventu- 
ally leading to a gross overestimation of the spin gap. 

In spin clusters E2-E5, all inter-dimer interactions are 
through the empty triangle. Similar to El, the ground 
state wave function is not given by a simple direct prod- 
uct of spin-singlet wave functions, but the overall spin- 
spin correlation is still well represented by the VBS order. 
To see this, look at the spin-spin correlation values given 
near each corresponding bond in Table |lTl Along the 
thick lines representing the valence bonds, the spin-spin 
correlation values are close to —0.75 J while their mag- 
nitude is much reduced along the bonds in the empty 
triangle structure. The ground state energy of empty 
triangle spin clusters is much lower than that of the topo- 
logically orthogonal counterparts, —0.3 75 J. In fact, the 
lowest energy configuration is obtained in E3 where three 
empty triangles form a closed circuit resulting in the per- 
fect hexagon. As for the excitation spectrum, low-energy 
spin excitations are quite sparse compared to those of 
the topologically orthogonal spin clusters. This result 
confirms again that the main role of the empty triangle 
structure is to lower the ground state energy. 

All spin clusters considered so far are entirely com- 
posed of the inter-dimer interaction structures in a single 
type. In Table IIIIl we investigate what happens when 
both types coexist. The simplest inter-dimer structure 
containing both structures, TE, is composed of a sin- 
gle topologically orthogonal structure, Tl, and a single 
empty triangle structure. El, as depicted in Table Hill 

One of the most important features of TE is that the 
ground state overall maintains the VBS order of each in- 
dividual inter-dimer structure. That is to say, being the 
highest in magnitude, the spin-spin correlation value is 
quite close to -0.75 along (1,2), (3,4), (5,6), and (7,8) 
link. Note that stands for the link connecting spin 
i and j. Meanwhile, the corresponding spin-spin corre- 
lation value is very small for the (5,7) and (6,7) link in 
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-0.426 


0.574 


10 






— (E3) 


-0.443 


0.542 


10 






(E4) 


-0.430 


0.518 


13 






(E5) 


-0.432 


0.487 


16 






(TE) 


-0.408 


0.392 


11 






^ In spin clusters studied, the lowest spin excitations are all 



(T1 ) (T2) (T3) (T4) (T5) (T6) ( E1 ) ( E2) ( E3) ( E4) ( E5) (TE) spin triplets. 

TABLE II. (Color online) Schematic diagram for the spin clusters containing only the topologically orthogonal inter-dimer 
interactions, T2-T6, and the empty triangle interactions, E2-E5. Dimers (light blue) indicate the valence bonds formed in the 
ground state. In the spin clusters T2-T6 where entire dimer connections are topologically orthogonal, the ground state wave 
function is precisely given by a direct product of the spin-singlet wave functions formed in each dimer. In this situation, the 
spin-spin correlation value, (Si • Sj), becomes exactly -0.75 within dimers and between dimers. It is important to note that 
each dimer configuration in T2-T6 is one of the twofold degenerate ground states. The other state not shown in the figure 
is obtained by reversing the orientation of the dimer covering pattern in the pinwheel structure. In the spin clusters E2-E5 
where all inter-dimer interactions are through empty triangle structures, the spin-spin correlation values (shown near each 
corresponding bond in the figure) of the ground state are overall similar to those of the simple spin-singlet product state while 
the wave function is not precisely in the product form. Lower left figure: energy spectrum of each spin cluster up to IJ of 
the excitation energy. In this figure, TE is a new spin cluster containing the both topologically orthogonal and empty triangle 
structure. See Table [Till for detail. Lower right table: the ground state energy per spin, tgr/J, and the spin-1 excitation 
energy gap, As=i/J. Nij denotes the number of states within IJ of the excitation energy from the ground state. 

12 3 4 



8 7 

(TE) 





El J 


5" 


(Si -Sa) 


(S3 ■ S4) 


(S5 ■ Se) 


(Sr ■ Sg) 


(Sa -Ss) 


(S3 ■ Ss) 


(Ss -Sa) 


(Ss - Sr) 


(Se • Sr) 


Ground state 


-3.265 





-0.676 


-0.676 


-0.662 


-0.742 


-0.134 


-0.169 


-0.169 


0.013 


-0.050 


1^* excited state 


-2.873 


1 


-0.623 


-0.623 


0.139 


-0.470 


-0.037 


-0.280 


-0.280 


-0.254 


-0.446 


2"'' excited state 


-2.535 


1 


-0.216 


-0.216 


-0.633 


-0.729 


-0.448 


-0.111 


-0.111 


-0.007 


-0.064 



TABLE III. (Color online) Energy, Ej J , Spin, 5", and the spin-spin correlation value, (Si ■ Sj), for the ground, the first-excited 
and the second-excited state in a spin cluster, TE, containing the both topologically orthogonal and empty triangle structure. 
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a close agreement with a single topologically orthogonal 
unit structure, Tl. Also, along the (2,3), (3,5), and (2,5) 
link, the spin-spin correlation value is roughly the same as 
that of a single empty triangle unit structure. El, —0.15. 
Another important feature of TE is that its topologically 
orthogonal part, Tl, has very little contributions to the 
lowering of the ground state energy. With the ground 
state energy per spin being —0.408 J, the contribution of 
Tl alone, J (S5 • Sg + S7 • Sg + S5 • S7 + Se • S7), is only 
—0.360 J per spin. It is due to the empty triangle struc- 
ture that the ground state has such a low energy men- 
tioned in the above. 

As for the excitations, it is found that the first excited 
state of TE has strong characteristics of the lowest-energy 
excitation in Tl while the second excited state is rather 
close to that of El. To be specific, the first excited state 
of TE is essentially obtained by breaking the spin singlet 
bond formed in the (5,6) link into a spin triplet, which 
can be seen via the change in the spin-spin correlation 
value from —0.662 to 0.139. In a close analogy with the 
lowest-energy excitation in Tl, this spin-singlet breaking 
in turn spreads the spin-singiet characteristics along the 
(7,8) link to nearby links of (5,7) and (6,7). Similarly, the 
second excited state of TE resembles the lowest-energy 
excitation of El. Note that this is consistent with the 
fact that the empty triangle structure has a higher spin 
excitation energy gap than the topologically orthogonal 
counterpart. 

To summarize, three lessens are obtained from exact 
diagonalization of various small spin clusters: (i) there is 
absolutely no quantum fluctuation about the VBS order 
in the spin clusters entirely composed of the topologically 
orthogonal structures, (ii) the empty triangle structure 
plays an important role for the lowering of the ground 
state energy, and (iii) the VBS order is quite robust in the 
ground state even when two different types of the inter- 
dimer structures coexist. This observation provides a mo- 
tivation for the VBS order in the Kagome lattice antifer- 
romagnet. In our picture, the strong geometric frustra- 
tion in the Kagome lattice prefers the ground state mini- 
mizing the quantum fluctuation. The VBS order achieves 
this goal by maximizing the number of topologically or- 
thogonal inter-dimer structures with a 36-site unit cell 
within which quantum fiuctuations are completely sup- 
pressed. This suppression, however, cannot occur in the 
whole lattice (unless the lattice is the Shastry-Sutherland 
lattice). Across the boundary of the unit cell, inter-dimer 
interactions are in fact solely mediated by the empty tri- 
angle structure. In this situation, the 36-site unit cell 
should be arranged in such a way that the ground state 
energy is minimized. As mentioned previously, the low- 
est energy configuration is obtained when three empty 
triangles are bound into a closed packet, leading to the 
perfect hexagonal structure. Two such final VBS order 
patterns are considered in this paper. 

In Sec. Illll and llVl we present a theoretical framework 
providing a natural platform for describing the VBS order 
ground state as well as its excitations. Before we do that. 



it is instructive to test how well the bond operator theory 
works for the two inter-dimer interaction unit structures, 
Tl and El. This test sets a stage for the spin cluster 
operator theory where energy eigenstates emerging from 
extended spin clusters, rather than those emerging from 
a single bond, are treated as independent quasi-particles 
of the system. 



B. Testing the bond operator theory 

We begin by writing the Heisenberg-model Hamilto- 
nian for a single topologically orthogonal unit structure, 
Tl, shown in Fig. [2J 

Ht1,J = J (S, • Sfc + Se • Sd) + A J {Sa + Sb) • Se, (2) 

where A is introduced for computational convenience, but 
is set to be unity at the end of computation. In the bond 
operator theory;^ the Hamiltonian is rewritten in terms 
of a new class of operators for collective quasi-particles 
emerging from a bond. 

Specifically, the bond operators are defined as follows: 





- |0) : 




= 4lo) 


ty) 


= 4io) 




= A |o) 



1 

71 



V2 

1 

71 



(in) + , 
(in) + in)) , 



(3a) 
(3b) 
(3c) 
(3d) 



where and ij, {a — x, y, z) are the respective operators 
for creating a spin singlet and triplet state. While either 
boson or fermion statistics can be assigned to bond par- 
ticles as long as the Lie algebra for the physical spin is 
correctly reproduced, it is convenient to treat both s and 
to, as boson operators. 

The spin operators for the left and right site, Sl and 
Sfl, are represented in the bond operator formalism as 
follows: 

SLa = ^ {sha + 4s - ieapff^i^tj^ , (4a) 
^(^-sha-As-ieaPjt^ptj^ , (4b) 



S 



Ra 



where h — 1, a, (3,^ G {x, y, z}, and the Einstein summa- 
tion convention is used. While correct, the above repre- 
sentation alone is not yet sufficient to define the physical 
spin operator. What is necessary is the hard-core con- 
straint: 



s - 



1 



(5) 



In other words, the spin operator representation in the 
above can generate the correct Lie algebra only in the 
Hilbert space restricted by the hard-core constraint. 
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The main reason why the bond operator theory is use- 
ful is that it provides a convenient framework to describe 
the VBS order. The VBS order is represented in the bond 
operator theory through the condensation of s bosons. 
TechnicaUy speaking, this means that both the annihila- 
tion and creation operator, s and s^, are replaced by a 
simple number, s, which denotes the condensation den- 
sity of spin-singlet pairs. 

With everything put together, the final Hamiltonian in 
the bond operator representation is written as follows: 



Hti — Hti,.j{s, to) 



2 

E 

i=l 



1 



(6) 



where i — \ and 2 distinguish between the two valence 
bonds in Tl. Here, //^ is the Lagrange multiplier used for 
imposing the hard-core constraint in the i-th. bond. Note 
that HTi,j{s,ta) is the Heisenberg-model Hamiltonian 
represented in terms of the s and ta bosons. Explicitly, 
the above Hamiltonian can be written as follows: 



Hti — £ti,o + Hti.2 + H^i^i, 



where 



--Js2+^,(l-S?) 



H' 



Tl,2 



E 

1=1 



t 



XJ 



t t 

"t1,4 = :T-£al3i£afj.utij3tl.yt2^t2,^. 



(7) 

(8) 

(9) 
(10) 



Note that all cubic terms containing three t operators 
are ignored in the above since they play little role in a 
non-magnetic ground state. It is interesting to observe 
that the quadratic Hamiltonian, Hti. 2, is actually inde- 
pendent of the inter-dimer interaction strength, A J. This 
means that the ground state is a simple direct product 
of two spin singlet states. As mentioned previously, this 
is consistent with the exact result for the ground state. 

The mean-field analysis is performed via quadratic 
decoupling of the quartic Hamiltonian, Htia- With 
two mean-field parameters defined by P = ('^la^Ja) ^^'^ 
Q — (iiQ^2a), the quadratically decoupled Hamiltonian 
of Hti,4 is given as follows: 

^quad. doc. ^ hl(Q*Q-P*p) 



3 
AJ 

~3~ 



{Pt\j2a~Qt\jL+ll. C.).(ll) 



For details of the derivation, see Appendix [Xj Putting 
back Eq. (|lip into (O , we obtain the following mean-field 
Hamiltonian: 



H^^ ^eri.o + ^{Q*Q-P*P) 
+ itlMt„ - ^TrM, 



(12) 



where 



M = 



f - Ml 

A J JZ>* J_ 



3 ^ 
- M2 















4 - Ml 

AJ p 

3 ^ 



M. P 

3 ^ 

J 



M2 / 



and 



/iia \ 

t2a 

V4J 



(13) 



(14) 



The mean- field Hamiltonian, H^^ , is diagonalized via 
the Bogoliubov transformation of ta ■ 



Hti = egr,Ti + E'^*^i«7ia, 



(15) 



where 



A 7 3 3 

egr,Ti = eTi,o + ^(Q*Q - P*P) - tTiM + ^ ^ u;. 



i=i 



(16) 

and u)i is the eigenfrequency of the i-th Bogoliubov quasi- 
particle, 7^^, which is in turn defined by 



J2a 
lie 



(17) 



with T being the Bogoliubov transformation matrix. 
Note that the Bogoliubov transformation in the above 
is different from that of the BCS theory where a fermion 
(electron) is transformed to another fermion quasiparti- 
cle. Here, a boson, ta, is transformed to another boson 
quasiparticle, 7q. To satisfy the canonical boson com- 



mutation relations, 



r r 



= lB<5a^, where Is is the 
diagonal matrix with elements being (1,1,-1,-1), the Bo- 
goliubov transformation matrix, T, is required to satisfy 
following constraint: 



TleTt = I, 



(18) 



Due to this constraint, the diagonalization of the boson 
problem is slightly different from the fermion counter- 
part. The boson eigenvalue problent^^ is given by 



TIbMT 



(19) 



where r2 is the diagonal matrix with elements being 

(wi,a;2,wi,W2)- 

Assuming that every Wi is non-negative for the stable 
valence bond ground state, |gr)rp^, the ground state en- 
ergy is equal to Cgr.Ti which is a function of A, s^, /x^, 
P, and Q. For a given A, Si and ^.i are self-consistently 
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determined through the fohowing four saddle point equa- 
tions together with P = {tiat\^) and Q = {tiat2a)'- 



gr.Tl 



dSi 
5egr,Tl 



4 

{1-s. 



J - fii] 2s,, = 0, 



(20a) 
0, (20b) 



where i = 1,2. The solution at A = is simple and given 

by 



1, 



Ml 



-J, P^Q, Q 0. 



(21) 



An important fact about the bond operator mean field 
theory is that the solution at A = remains to be the 
solution at any A up to unity. This means that the ground 
state is simply given by the direct product of two spin 
singlet states: 



|gr; 



Tl 



/ab ' 



I cd ' 



(22) 



where |s)j^- denotes the spin singlet wave function formed 
between the i and j-th site. Excited states are induced by 
breaking either the (a, b) or the (c, d) valence bond into 
a spin triplet, which results in a degenerate excitation 
energy of IJ. 

It is important to note that, in comparison with the 
exact states, the first excited state turns out to be grossly 
inaccurate while the ground state is correctly captured by 
the bond operator mean field theory. See Table IIVI for 
details of the result. In the bond operator mean field 
theory, the wave function for the first excited state is 
given by 



/ cd ' 



(23) 



7Llgr)Ti = K)ab' 
while the exact solution at A = 1 is 

^ iK)ab ® k>cd + \to.)bd ® l«>ca + \tc.)ad 1^) cb) ■ 

(24) 

As one can see, similar to the bond operator result, the 
exact solution for the first excited state has a broken va- 
lence bond in the (a, b) link. However, due to additional 
spin correlations emerging from the (a, d) as well as (fo, d) 
link, the final spin-spin correlation appears so that the 
spin singlet, which is supposed to form in the (c, d) link, 
spreads into the (a,c), (6,c), and (c, d) link equally with 
the spin-spin correlation value becoming —0.417. Conse- 
quently, the excitation energy of the exact first excited 
state, 0.5 J, is much lower than that of the bond operator 
mean field theory, IJ. Interestingly, the second excited 
state is well captured by the bond operator mean field 
theory. To summarize, while correctly representing the 
ground (and the second excite) state, the bond operator 
theory description is rather inaccurate for the lowest spin 
excitation in the topologically orthogonal unit structure, 
Tl. The conclusion so far is that the simple bond oper- 
ator theory is not accurate enough to correctly capture 



State 


Spin 


Energy 


{a,b) 


(a,c), (6,c) 


(c,d) 


!gr)Ti 





-1.5J 


-0.75 





-0.75 




(—0 75) 


fo) 


f-O 75) 


7L lgr>Ti 


1 


-0.5J 
(-1.0J) 


0.25 
(0.25) 




(-0.417) 


-0.75 
(-0.417) 


7L |gr>Ti 


1 


-0.5J 
(-0.5J) 


-0.75 
(-0.75) 




(0) 


0.25 
(0.25) 



TABLE IV. Various properties of the low energy states in the 
topologically orthogonal unit structure, Tl, obtained from the 
bond operator mean field theory. The exact diagonalization 



result is shown in 


parenthesis 


for comparison. 


State Spin 


Energy 


(a, b), (c, d), (e, /) (6, d), (d, e), (e, b) 


|gT)El 


-2.524J 
(-2.5J) 


-0.666 -0.176 
(-0.683) (-0.15) 


7mcK |g^)El 


-1.629J 
(-1.764J) 


-0.572 0.029 
(-0.370) (-0.218) 



TABLE V. Various properties of the low energy states in the 
empty triangle unit structure, El, obtained from the bond 
operator mean field theory. The exact diagonalization results 
are shown in parenthesis for comparison. Note that the first 
excited states are sixfold degenerate and distinguished by m = 
— 1,1 and a = x,y, z. 



the lowest spin excitation. We believe that a promis- 
ing direction for improvement is to construct a theory 
that incorporates the spin-spin correlations of the entire 
topologically orthogonal structure, not just between two 
nearest spins inside. The spin cluster operator theory 
(SCOT) is such an attempt. 

Now, let us move to the empty triangle unit structure. 
El, whose schematic diagram is shown in Fig.[2j Relegat- 
ing calculation details to Appendix [B] here, we discuss 
the results of the bond operator mean field theory and 
compare them with exact diagonalization. First of all, as 
one can see in Table |Vl the ground state energy as well as 
the spin-spin correlation values are in an excellent agree- 
ment with the exact diagonalization results. This means 
that the ground state is quite well described by the bond 
operator mean field theory even though it is not the exact 
eigenstate. 

As for the lowest excitation, the bond operator the- 
ory predicts sixfold degenerate spin excitations, which is 
consistent with the exact results. Moreover, while there 
are some discrepancies between the spin-spin correlation 
value of the bond operator mean field theory and exact 
diagonalization, the energy itself is very well reproduced 
by the bond operator mean field theory. Considering 
the excellent agreement in energy, we believe that the 
bond operator mean field theory works reasonably well 
in the case of the empty triangle structure. It is impor- 
tant to note that the inter-dimer interaction is completely 
ignored by the bond operator mean field theory in the 
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topologically orthogonal structure while it is not true in 
the empty triangle counterpart. 

The final conclusion of this section is that the ground 
state is quite well captured by the bond operator mean 
field theory in the both topologically orthogonal and 
empty triangle structure. The lowest energy excitation 
is, however, poorly represented in the topologically or- 
thogonal structure since the simple bond operator the- 
ory regards a pair of the nearest neighboring spins as the 
basic building block. We believe that it is necessary to 
extend the basic building block of the system dynamics 
from a dimer to a spin cluster to incorporate as many 
topologically orthogonal structures as possible. 



C. Binding dimers into a spin cluster 

As mentioned in the preceding section, for the correct 
representation of low-energy spin excitations, it is neces- 
sary to bind dimers into a spin cluster so that the spin- 
spin correlation in the topologically orthogonal structure 
is better taken into account. To distinguish from the con- 
ventional bond operator theory, we call such theory the 
spin cluster operator theory (SCOT). 

As depicted in Fig. [31 there are two possible patterns 
for the VBS order with a 36-site unit cell: the VBS order 
with (i) staggered hexagonal resonance and (ii) uniform 
hexagonal resonance. Appropriate spin clusters are se- 
lected according to the following conditions: (i) the VBS 
order formed inside each spin cluster should be consistent 
with that of the whole lattice, (ii) the shape of the spin 
clusters should preserve the symmetry of the VBS order, 
and (iii) the spin clusters should include as many topo- 
logically orthogonal structures as possible within them- 
selves. The most general spin cluster that satisfies all 
these three conditions is in fact the unit cell itself, where 
every inter-dimer interaction is topologically orthogonal. 
Dubbed as the "super cluster," however, the whole unit 
cell is too large to be diagonalized directly since a very 
large number of spin triplet excited states are necessary 
to correctly represent the physical spin. In this study, we 
choose a particular subset of the unit cell that is consis- 
tent with all the above three conditions. 

Figure [3] shows the spin cluster covering for the VBS 
order with staggered and uniform hexagonal resonance. 
Let us first consider the VBS order with staggered hexag- 
onal resonance, in which case the unit cell is composed 
of three spin clusters: a David star, D, and two extended 
hexagons, HI and H2. The David-star spin cluster con- 
sists of six dimers which form a pinwheel structure. Ex- 
tended hexagons are created by binding the rest of the 
dimers in such a way that the resulting spin cluster cov- 
ering preserves the sixfold (threefold) rotation symme- 
try around the center of the pinwheel (perfect hexagon) 
structure. It can be shown that the above separation of 
the unit cell into subset spin clusters is fully consistent 
with the VBS order of the whole lattice. For the VBS or- 
der with uniform hexagonal resonance, there are two dif- 



-0.5 

7 
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FIG. 4. (Color online) Schematic diagram for the David-star 
spin cluster and the excitation energy spectrum showing two 
thousands states above two degenerate ground states at S = 
0. 

ferent kinds of the spin clusters in the unit cell: (i) a core 
cluster, C, composed of a David star and three attached 
dimers and (ii) nine dimers, SI - S9. It is confirmed again 
that such separation into subset spin clusters is consistent 
with the VBS order with uniform hexagonal resonance, 
which has a threefold rotation symmetry around the cen- 
ter of the pinwheel structure and a reflection symmetry 
with respect to the vertical axis piercing through the pin- 
wheel structure. Now that appropriate spin clusters are 
selected for each VBS order, we develop a fully hedged 
spin cluster operator theory in the next section. 



III. SPIN CLUSTER OPERATOR THEORY: 
STAGGERED HEXAGONAL RESONANCE 

A. Spin cluster operator representation 

1. David star 

Constructing the spin cluster operator representation 
is actually nothing but choosing an appropriate basis set 
suitable for representing the low energy states in the pres- 
ence of the VBS order. In choosing a suitable basis set 
for the spin cluster, first noted is the fact that the num- 
ber of available states increases exponentially as the size 
of the spin cluster grows. For example, in a David star 
spin cluster, there are 2^^ states with various spin values 
ranging from = to 6. The figure |4] shows two thou- 
sand low energy states emerging from a David star spin 
cluster. As a matter of principle, all eigenstates should 
be retained for the exact representation of the physical 
spin. Setting aside the technical issue that the number of 
eigenstates is too large, an important physical question 
is how to deal with the large spin states. In this work, we 
take an approach that, first, the lowest energy state in 
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the S = channel is condensed so that the appropriate 
VBS order is formed in the ground state. Second, with no 
clear energy gap in the S = 1 channel dividing between 
the low and high region of the spectrum, we include as 
many spin triplet states as possible for the representation 
of the physical spin. Technically, our approach is to in- 
crease the number of spin triplet states until results such 
as the ground state energy as well as the spin gap con- 
verge into a well-defined fixed value. Finally, assuming 
the the lowest spin excitations are spin triplet, i. e., in 
the S = 1 channel, we ignore all other higher spin states 
for the spin cluster representation. 

As mentioned in the preceding section, the ground 
state in the spin singlet channel is degenerate with two 
different orientations for the dimer covering pattern in 
the pinwheel structure. Mathematically, the ground- 
state Hilbert space is spanned by the following two VBS 
states: 



101 



102) = 



(25a) 



(25b) 



(«j)e/2 



where /i={(l,8), (2,9), (3,10), (4,11), (5,12), (6,7)} and 
/2={(1,7), (2,8), (3,9), (4,10), (5,11), (6,12)}. The above 
two states are not orthogonal as seen from the non-zero 
overlap: (0i|(/)2) — 1/32. While it is straightforward to 
construct an orthonormal basis set out of the above two 
states: 

IV-i) = ^=i==(|0i) + |02)), (26a) 



\/2 + 1/16 



1^-2 



^2-1/16 



(I0l)"l02)), 



(26b) 



the new states, jV'i) and 1-02), no longer possess the VBS 
order. In light of our assumption that the VBS order 
exists in the ground state, we take either |0i) or |02) to 
be condensed. It does not matter which one of the two 
states is condensed since the reflection symmetry con- 
necting between |0i) and \<p2) is presumed to be broken 
spontaneously. 

As for the excited states in the S' = 1 channel, ba- 
sis states are chosen first as the eigenstates of (S = 
J2i=i^i and Sa (ck = x,y,z) in the following combina- 
tion: 

\t^) = ^ (|5 = 1, TO = 1) - |5 = 1, m = -1)) , (27a) 



\ty) = — i\S=l,m = l) 

\t,) = 1,TO = 0). 



1, m 



-1)), (27b) 
(27c) 



Note that each \ta) is invariant under their respective 
spin rotation, i. e.. 



Sa; \tx) — Sj^ \ty) 



0. 



(28) 



Unfortunately, the above classification does not uniquely 
choose the basis states since there exist degeneracies 
among many of the spin triplet states. This problem 
can be resolved by considering the sixfold space rotation 
operator, Rg, which is defined as follows; 

Re |toi, TO2, 1713,1714, mz, me, tuj, ms, TOg, mio, mn, mi2) 

= |to6, mi, TO2, ms, TO4, ms, toi2, tot, ms, toq, mio, tou) , 

(29) 

where rrii denotes the eigenvalue of Si,z. The eigen- 
value of Rg is given by an element of the set, Mg = 
|gi27r;/6|; ^ -2,-1,0,1,2,3}. Degenerate spin triplet 
states can be decomposed into the eigenstates of Rg since 
it commutes with both the Hamiltonian, ' 
and the square of the total spin operator, S^. Degener- 
acy in energy can appear between Vs with different signs, 
but the same magnitude, i. e., either / = ±1 or I = ±2. 

Taking the eigenstates of Rg provides an additional 
advantage for the construction of the spin cluster repre- 
sentation. If it were not for Rg, the spin cluster repre- 
sentation has to be computed for every individual spin. 
With help of the sixfold space rotation symmetry, all spin 
operators can be computed from the knowledge of only 
two spin operators, Si and S7: 



Si+p — RgSiRg 



(30a) 
(30b) 



for p = 0, . . . , 5. 

The spin cluster operator representation of the physical 
spin is formally written as follows: 



Si+p,Q — ^^Ol (^1 RgSiaRg 



(31) 



where i — 1,7 and p = 0, • • • , 5. Here, is the fi- 
th energy eigenstate and denotes the corresponding 
creation operator. Restricted to the Hilbert space with 
the ground state (which has a VBS order) and the spin 
triplet excitations, the spin cluster operator representa- 
tion of the physical spin is finally obtained as follows: 



S, 



+p,a ■ 
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m,n— 1 
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n=l 
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mn 



(32) 



where is the creation operator for the spin singlet 
ground state with one of the two possible pinwheel con- 
figurations, |0i) and \(t>2)i in Eq. 4ia is the creation 
operator for the spin triplet state which is a simultane- 
ous eigenstate of the Hamiltonian, H , the sixfold space 
rotation operator, Rg, and the spin operator, S^: 

H \tna) — l^ria) 7 (33a) 
R.6 \tna) = Zn \tna) {Zn G ^g), (33b) 
Sa |t„a) = 0, (33c) 
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and Sg, determine those of other eight spins with help of 
the threefold rotation symmetry: 



+2p 



(37) 



where i = 1, 2, 7, 8 and p = 1,2. 

Repeating essentially the same procedure described in 
the previous section, the spin cluster operator represen- 
tation for the extended hexagon is obtained as follows: 



-*i-\-2p,a 
N 



N 

E 

n=l 



H. c. 



(38) 



FIG. 5. (Color online) Schematic diagram for the extended 
hexagon spin cluster and the excitation energy spectrum 
showing two thousands states above the spin-singlet ground 
state. 



where is the eigenenergy of the n-th spin triplet ex- 
citation in the David star and z„ is the eigenvalue of Rg 
belonging to Mg. In the above, M is the total number 
of spin triplet states included in spin cluster operator 
representation. The matrix elements, fn 
defined by 



where 



and e 



{^my I S 



(i) 



(34) 
(35) 



where it is interesting to note that e,n„ is anti-hermitian. 
A brief sketch for the derivation of the spin cluster oper- 
ator representation is given in Appendix [Cj 



2. Extended hexagon 

As mentioned previously, the unit cell of the VBS 
ground state with staggered hexagonal resonance is com- 
posed of a David star and two extended hexagons with 
three attached dimers. The spin cluster operator repre- 
sentation for the extended hexagon is constructed in the 
similar way to the David star case with two minor mod- 
ifications. First, the ground state has no degeneracy in 
the extended hexagon spin cluster. Second, the extended 
hexagon spin cluster has a threefold space rotation sym- 
metry. With the threefold space rotation operator, R3, 
defined by 

R3 |toi, m2, TOa, 7714, "T-s, "76, mj, TTjg, 7779, mio, mn, 77712) 

= I7775, 7716, mi, 7772, r773 , 77l4, 777ii , 777 i2 , 77l7, 7778, 7779, mio) , 

(36) 

the eigenvalue of R3 is an element of the set, K3 = 
|gi2ii-V3|^ _ _i^o, 1}. As before, the spin cluster opera- 
tor representation of the four spin operators. Si, S2, S7, 



R3 \tna) 



rn \tna) (^i € M3), 
(^nx I ^ix I ^} 1 
{^7ny I ^ix l^nz) 5 



(39) 
(40) 

(41) 

(42) 



and N is the total number of spin triplet states considered 
for the extended hexagon spin cluster. 



B. Hamiltonian and the mean field theory 

The basic plan for analysis is to rewrite the Heisenberg- 
model Hamiltonian in the spin cluster operator represen- 
tation, which is obtained in the preceding section, and 
apply the self-consistent mean-field method similar to the 
bond operator theory as described in Section lll Bl Let us 
begin by formally writing the Heisenberg-model Hamilto- 
nian in the spin cluster operator representation, accom- 
panied by three Lagrange multiplier terms each corre- 
sponding to the hard-core constraint in a David star and 
two extended hexagons: 

M 



H = HJ-^,nJ2T. [^D + Cnr)C(r) - l] 

r n— 1 

TV 

-/^HEE[«"l^+^«i'«C^w-i] 

r n—1 

N 

- E E I'h + (r) - 1] , (43) 

r n—1 

where 

Hj = J S, -Sj+AJ E S,-S,, (44) 
{i.j)ec {i.j)^c 

with C denoting the site-index pairs belonging to the va- 
lence bond in the David star, D, as well as the two ex- 
tended hexagon spin clusters, HI and H2. As before, 
A is introduced for computational convenience and is set 
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to be unity at the end of computation, s and /i denotes 
the spin singlet condensation density and the Lagrange 
multipher for the hard-core constraint, respectively. 

Arranging the Hamiltonian with respect to the degree 
of spin triplet operators leads to 

H = ^eo + H2+ Ha, (45) 

where ,yV is the total number of unit cells. H2 and 
is the quadratic and quartic part of the Hamiltonian, 
respectively. The constant term in the Hamiltonian, eo, 
is given by 



'^0 



(46) 



where €q and is the ground state energy of the David 
star and the extended hexagon spin clusters, respec- 
tively. The mean field Hamiltonian is obtained via self- 
consistent quadratic decoupling of the quartic terms as 
described in Appendix |^ Mean- field parameters are de- 
fined as follows: 



kl 



^DH2 
kl 



(r') 
(r') 



kl 



(C(r)tr(r + r') 



C(r)C^(r + r') 
C(r)C^(r + r'; 



and 



^r^(r') 

^.r(r') 



(C(r)C(r 
(C(r)^S^(r 
(C(r)C(r + r')), 
with r' being one of the following vectors: 




4%/3a (-y) 



(47a) 
(47b) 
(47c) 

(48a) 
(48b) 
(48c) 

(49a) 

(49b) 
(49c) 



where a is the distance between the nearest neighboring 
spins. 

The mean field Hamiltonian can be written in a com- 
pact form via the Fourier transformation: 



MF 



2 ^ 

k 



Ti(k)M(k)r„(k)~^VTrM(k), 



4 ^ 

k 



(50) 



where 



Ta(k) = 



/ t^(k) \ 

t^^(k) 
tg^(k) 

en-k) 

Vtf^t(-k)y 



(51) 



In the above, tj^(k) and tQ(k) are the respective column 
vectors composed of the creation and annihilation oper- 
ators for spin triplet states. The size of the r vector 
is given by 2(M -I- 2A^). The detailed form of and 
M(k) are given in Appendix ID] 

Diagonalizing the mean field Hamiltonian via the Bo- 
goliubov transformation leads to the following ground 
state energy per spin: 

1 , 

Cgr — ^ [ eo + e,^^ 

-!:^E™(k) + |^5: -"(k)], 



2 J/ 



k n=l 



(52) 



where w„(k) is the solution of the eigenvalue problem: 
T(k)lBM(k)T-i (k) = iBrj(k), (53) 



where 



Is = 
n(k) = 



I 
I 

a;(k) 
a;(k) 



(54) 
(55) 



As before, the spin singlet condensation density as well 
as the Lagrange multiplier for the hard-core constraint 
are determined by the saddle-point equations: 



gr 



ds 



0, 



(56) 



where s = sd or sh, and /i = /i/? or ^h- The saddle- 
point equations are solved simultaneously under the 
self-consistent conditions for the mean field parameters 
defined in Eq. (gT]) and (gSD. With a solution of the full 
parameter set, (s, /x, simultaneously satisfying 

all the saddle-point and the self-consistence equations 
at A = 1, one can compute the ground state energy as 
well as the spin triplet excitation spectrum. Results are 
discussed in Section IVl 



IV. SPIN CLUSTER OPERATOR THEORY: 
UNIFORM HEXAGONAL RESONANCE 

A. Spin cluster operator representation 

The core cluster shown in Fig. [B] has a twofold de- 
generacy in the ground state energy, which is related to 
two possible dimer orientations in the pinwheel structure. 
That is to say, either of the following two spin singlet 
states can be condensed for the ground state: 



(ij)e/uJi 



(57a) 
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FIG. 6. Schematic diagram for the core cluster in the VBS 
order with uniform hexagonal resonance. 

k2)= (8) (57b) 
(ij)e/uJ2 

where I={(9,1), (3,4), (6,7)}, Ji={(ll,2), (12,13), 
(14,5), (15,16), (17,8), (18,10)}, and J2={(10,11), (2,12), 
(13,14), (5,15), (16,17), (8,18)}. 

By using the similar method for the analysis of the 
David star as well as the extended hexagon cluster, the 
spin cluster operator representation for the core cluster 
is obtained as follows: 

L 
n=l 

L 

+ 5^ {^mn ' fm^n^)^al37tml3^ri-f, (58) 
m,n—l 

where i = 1, 2, 3 and p = 0,1, 2. Also, in the above, 

k)e{ki),|7r2)}, (59) 
R3 \s)^\s), (60) 

R3|tna)=r„|i„„)(f„eM3), (61) 

a^^ = {tn.\S^,\s), (62) 

^rnn ~ i^rnyl Six \tnz) j (63) 

and L is the total number of spin triplet states included 
for the core cluster representation. 

In addition to the core cluster, there remain nine 
dimers surrounding the core. See Fig. [3] for the illus- 
tration of the VBS unit cell with uniform hexagonal 
resonance. The nine surrounding dimers can be repre- 
sented by the usual bond operator formalism discussed 
in Sec. Hrei 

B. Hamiltonian and the mean field theory 

The analysis of the mean field Hamiltonian for the VBS 
order with uniform hexagonal resonance is similar to that 
of the staggered hexagonal resonance counterpart. Since 
computational details are basically identical, we skip the 
technical discussion and just mention a few points that 
are specific to the VBS order with uniform hexagonal 



resonance. First, once the core cluster is isolated from the 
unit cell, there is no natural spin cluster with an extended 
size compatible with the VBS order. For this reason, all 
remaining spins are paired as dimers and analyzed via 
the bond operator theory. Second, it turns out that the 
core cluster is completely decoupled from the surrounding 
dimers so that spin triplet excitations are separated into 
two independent channels with one being confined into 
the core cluster and the other propagating through the 
surrounding dimers. Details of the results are given in 
Sec. ED 



V. RESULTS: STAGGERED HEXAGONAL 
RESONANCE 

The results of the spin cluster operator theory for the 
VBS order with staggered hexagonal resonance are dis- 
cussed in this section. In addition to the thermodynamic 
limit where the site number goes to infinity, we also re- 
port the results for the 36-site unit cell, which can be 
directly compared with the exact diagonalization results 
obtained in the same-sized system. Here, the 36-site unit 
cell means that the periodic boundary condition is used 
without the appropriate phase factors required for the 
infinite lattice. 

In the both cases of the infinite lattice and the 36-site 
unit cell, computations are performed as a function of the 
number of spin triplet states included in the spin cluster 
operator representation. With the number of included 
spin triplet states denoted as M for the David star and N 
for each extended hexagon cluster, the total spin-triplet 
state number is increased along M = in this study. 
It is found in the end of computation that the ground 
state energy as well as the lowest spin excitation gap 
are quite well converged already for M ~ 80. Note that 
various parallel computation techniques are used in this 
study to cope with the complexity associated with a large 
number of self-consistent mean field parameters defined 
in Eq. (|47l) and (p ]) . 



A. Ground state energy 

The ground state energy per spin is shown as a func- 
tion of the spin-triplet state number included for the spin 
cluster operator representation in Fig. [7l which also con- 
tains a similar plot for the spin triplet excitation gap. 
Figure [7] shows computational results for four different 
situations: (i) quadratic only approximation ignoring all 
quartic terms in the 36-site unit cell, denoted as "36-site 
quad-only" in the figure, (ii) self-consistent mean field 
decoupling of the quartic terms in the 36-site unit cell, 
denoted as "36-site mean-field", (iii) quadratic only ap- 
proximation in the infinite lattice, denoted as "oo-lattice 
quad-only" , and (iv) self-consistent mean field decoupling 
in the infinite lattice, denoted as "oo-lattice mean-field" . 
As mentioned previously, the results begin to show a con- 
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M ( =N ) 

FIG. 7. (Color online) The ground state energy per spin, 
Egr/J, (top panel) and the lowest spin triplet excitation en- 
ergy gap, As=i/ J, (bottom panel) as a function of spin triplet 
states included for the spin cluster operator representation. In 
this figure, the number of the included spin triplet states for 
the David star, M, and for each of the extended hexagons, 
A'', are set to be equal and increased from 10 to 140. 



vergent behavior as the number of included spin triplet 
states in each spin cluster goes beyond about 80. 

Table I VII provides the fully self-consistent mean-field 
results at M — N = 140 for the 36-site unit cell, de- 
note as SCOT lag, and the infinite lattice, denoted as 
SCOT loo- Note that, while the 36-site cluster has a lit- 
tle bit lower ground state energy and a slightly higher 
spin triplet gap than the infinite lattice, the difference 
is very small compared to their absolute values. In the 
spin cluster operator theory, the ground state energy is 
estimated to be —0.428 J per spin for the 36-site unit cell 
while it is given as —0.438 J from exact diagonalization. 
Considering the crude nature of the mean-field approx- 
imation, we believe that such difference in the ground 
state energy is not too bad. Note that the energy differ- 
ence is only 2.3% of the ground state energy itself. 

Also, it is shown in Table IVll that the SCOT loo ground 
state energy is lower than that of the BOT loo by 0.006 J, 
which means that the spin cluster operator theory does a 
better job in utilizing the inter-dimer interaction to min- 
imize the energy. In fact, the BOT Iqo simply reduces 
to the bond operator theory for a single spin cluster, E3, 
described in Table HH This is because, in the bond op- 
erator mean field theory, any topologically orthogonal 
inter-dimer link behaves as if the link itself is absent. 



B. Spin triplet excitation spectrum 

Figure [8] shows the spin triplet excitation spectrum 
of the spin cluster operator theory for the VBS order 
with staggered hexagonal resonance. As mentioned pre- 
viously, the results presented in Fig. [8] are those obtained 
at the full convergence occurring when the total num- 
ber of the included spin triplet states for the spin clus- 
ter operator representation reaches M + 2A^ ~ 420 with 
M = N = 140. 

Before discussing details of the spin triplet excitation 
spectrum, it is convenient to introduce the following 
quantity characterizing the regional origin of each spin 
triplet excitation: 

(64) 

where I — 1, - ■ ■ ,M + 2N and T^^(k) is the inverse Bo- 
goliubov transformation matrix. Noting that the spin 
triplet excitations from the David star are enumerated 
with indices between 1 and M and those from the ex- 
tended hexagons are with indices between M + 1 and 
M + 2N, the above quantity measures which region, the 
David star or the extended hexagons, is more respon- 
sible for a given spin triplet excitation. In the above 
definition, anomalous components of the eigenvectors, 
[T-^k)U {n = {M + 2N) + \,-- - , 2{M + 27V)), are 
not included due to their little contribution. Each dis- 
persion curve in Fig. [5] is color-coded based on the value 
of C;(k) with the blue indicating that the correspond- 
ing spin triplet excitation comes from the David star and 
with the red otherwise. The yellow color denotes that 



TABLE VI. The ground state energy per spin, egr/ J, and the 
lowest spin triplet excitation energy gap, As=i/J, obtained 
from various computations. The SCOT I and II denote the 
results obtained from the spin cluster operator theory for the 
VBS order with staggered and uniform hexagonal resonance, 
respectively. Results from the 36-site unit cell and the infi- 
nite lattice are distinguished by the subscripts under I and II. 
Similarly, the BOT I and II address the bond operator the- 
ory for the corresponding VBS orders. Note that, while the 
BOT IIoo is essentially identical to the previous work done 
by some of the authorsi^,^ there is a slight technical difference 
in the quadratic-decoupling procedure resulting in a minor 
quantitative modification. 
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0.2431 


DMRG^ 




0.05 


BOT U 
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FIG. 8. (Color online) Spin triplet excitation spectrum obtained from the spin cluster operator theory for the VBS order with 
staggered hexagonal resonance. Figure (a) shows the triplet excitation spectrum with energy up to 1.5J. Each spin triplet 
dispersion is color-coded depending on which region between the David star and the surrounding extended hexagons is more 
responsible for the origin of the corresponding spin triplet excitation. The blue color indicates that the spin triplet excitation 
comes from the David star while the red color denotes otherwise. The yellow color indicates that there is an equal mixing 
between contributions from the both regions. Figure (b)-(d) show an enlarged view of the spectrum for the energy window 
of oj/J = 0.32 — 0.34, 0.8 — 0.9, and 1.3 — 1.4, respectively. Figure (e) shows the plotted path of the spin triplet excitation 
dispersion in the first Brillouin zone. 



two contributions are almost equal. 

As one can see from Fig.|8](a), most spin triplet excita- 
tion dispersions are nearly flat with few exceptions of the 
dispersive modes in high energy. While there is in general 
non-zero mixing between the spin triplet modes originat- 
ing from the David star and the surrounding extended 
hexagons, low energy spin excitations mostly come from 
the David star which is more or less disconnected from 
the outside as indicated by the blue color. In high energy, 
only a few dispersive modes, denoted in red, have dom- 
inant contributions from the extended hexagons. It is 
interesting to note that the lowest spin triplet dispersion 
in red are twofold degenerate and occurs around 0.75J 
with a flat dispersion. A similar behavior is observed in 
the case of the VBS order with uniform hexagonal reso- 
nance as shown in next section. 

In addition to the blue and red spin triplet excitations 
exclusively originating from either region of the David 
star or the extended hexagons, spin triplet excitations 
in yellow have the equal contribution from both regions, 
appearing between 0.8 and 0.9 J. Figure [5] (c) shows an 
enlarged view of the spin triplet excitation spectrum in 
such region. Roughly speaking, the spectrum becomes a 
superposition of the dispersive red and the relatively flat 
blue curve only when the two colored bands intersect each 
other. Away from the intersecting points, the spin triplet 
excitations possess mostly the single color characteristics. 

The lowest triplet band is contained in the Fig. [S] (b) 
together with the next four bands. The lowest spin triplet 
excitation occurs at the K point with energy of 0.323 J. 
It is interesting to compare the SCOT result for the low- 
est spin excitation with that of the usual bond operator 
theory. The spin cluster operator theory predicts a much 



smaller spin triplet excitation gap than the bond opera- 
tor theory. See Table IVII for details. Despite the better 
agreement, there still remains a factor of 3 — 6 difference 
in comparison with the numerical results obtained from 
exact diagonalization as well as the density matrix renor- 
malization group method. While this difference could be 
induced by the crude approximation in the mean field 
theory itself, we believe that it reveals a more fundamen- 
tal limitation of the current spin cluster operator theory. 
As mentioned in the above, the lowest spin excitation 
mostly originates from the David star. Expecting that 
the true low-energy spin excitation has more or less the 
same contributions from both the David star and the 
extended hexagons, we believe that it is necessary to de- 
vise a theoretical framework where the spin triplet exci- 
tations from the both contributions are treated on a more 
equal footing. We discuss such theoretical framework in 
SecEni 



VI. RESULTS: UNIFORM HEXAGONAL 
RESONANCE 

A. Ground state energy 

The ground state energy per spin for the VBS or- 
der with uniform hexagonal resonance is obtained as 
— 0.422J in the fully self-consistent mean-field approxi- 
mation, reduced from the quadratic-only approximation 
value, -0.414J. The SCOT result is in fact precisely 
identical to that of the bond operator theory as shown in 
Table IVTl The reason is that, in the VBS oder with uni- 
form hexagonal resonance, the core cluster is completely 
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theory, the actual lowest spin triplet excitation occurs 
from the core cluster instead of the surrounding dimers. 
For a similar reason mentioned previously in the case of 
the VBS order with staggered hexagonal resonance, we 
think that this separation between the core cluster and 
the surrounding dimers might be artificial. To overcome 
this problem, in the next section, we propose a new the- 
oretical framework for the spin cluster operator theory, 
where the both regions are treated in a more equal foot- 
ing. 



VII. DISCUSSION 



FIG. 9. (Color online) Spin triplet excitation spectrum ob- 
tained from the spin cluster operator theory for the VBS order 
with uniform hexagonal resonance. The spectrum consists of 
two separate sets of the spin triplet excitation with one com- 
ing from the core cluster and the other from the surrounding 
dimers. The second lowest flat bands from the surrounding 
dimers are twofold degenerate. It is interesting to note that 
similar flat bands appear in the VBS order with staggered 
hexagon resonance with roughly the same energy. 



decoupled from the surrounding dimers. It is mentioned 
in Sec. Ill Al that such decoupled core cluster does not have 
any quantum fluctuation about the pure dimer-covering 
VBS order and therefore the lowest energy per spin for 
the core cluster is simply given by the half energy of a spin 
singlet, —0.375 J. Remaining contributions to the ground 
state energy are generated through the inter-dimer inter- 
action between the surrounding dimers. These contribu- 
tions, however, are exactly the same between the spin 
cluster and the bond operator theory. 



B. Spin triplet excitation spectrum 

The spin triplet excitation spectrum for the VBS order 
with uniform hexagonal resonance consists of two sepa- 
rate sets of excitations with one coming from the core 
cluster and the other from the surrounding dimers. See 
Fig. [S] for the spectrum of each set. Spin triplet excita- 
tions coming from the core cluster are entirely flat in dis- 
persion due to the fact that the core cluster is completely 
decoupled from the surrounding dimers. It is shown in 
Fig. m that the lowest spin triplet excitation has an exci- 
tation energy gap of 0.243 J. 

The surrounding dimers provide nine spin triplet ex- 
citations with five of them being fiat in dispersion. Be- 
longing to the five fiat modes, the lowest spin triplet ex- 
citation has an energy of 0.776 J and is degenerate at the 
F point with a dispersive mode. Similar to the ground 
state energy, these results are identical to those of the 
bond operator theory since the connection between the 
core cluster and the surrounding dimers is absent. The 
difference is that, in the current spin cluster operator 



The spin cluster operator theory predicts a spin triplet 
excitation gap significantly reduced from that of the bond 
operator theory. This reduction results from the fact that 
effects of the topologically orthogonal inter-dimer inter- 
action, which were entirely ignored in the bond operator 
theory, are better captured by the spin cluster operator 
theory. Despite the improvement, however, the current 
spin cluster operator theory is not accurate enough to 
correctly describe the full effects of the topologically or- 
thogonal inter-dimer interaction between dimers across 
different spin clusters. One way to see the inaccuracy is 
to realize that the spin triplet excitation gap is actually 
increased after the inter-dimer interactions between the 
David star and the extended hexagons are included in the 
VBS order with staggered hexagonal resonance. This is 
in contradiction with the observation that spin clusters 
with more topologically orthogonal inter-dimer interac- 
tions generate a smaller spin triplet excitation gap. 

An apparent source of the problem is that spin clusters 
considered in our current formalism do not fully take into 
account the effects of all topologically orthogonal inter- 
dimer interactions. A natural conclusion based on this 
reasoning is that the best way to overcome this problem 
is to construct a "super cluster" that contains all topolog- 
ically orthogonal inter-dimer interactions inside the unit 
cell. See Fig. [TU] for the illustration of possible super 
clusters. Note that, similar to the David star in the VBS 
order with staggered hexagonal resonance and the core 
cluster in the uniform hexagonal resonance, two super 
clusters corresponding to each hexagonal resonance have 
twofold degeneracies with respect to the dimer covering 
orientation in the central pinwheel structure. 

It is mentioned in Sec. Ill Al that there is absolutely 
no quantum fluctuation in the VBS state when all inter- 
dimer interactions are topologically orthogonal. That is, 
the VBS state is a precise energy eigenstate of a spin 
cluster entirely composed of the topologically orthogonal 
inter-dimer interactions, in which case the energy of the 
VBS state is exactly given by —0.375 J per spin. Actually, 
it is shown in Table [Tl] that the VBS state is not only a 
precise energy eigenstate, but also the exact ground state 
for all finite-size spin clusters studied in this paper, (Tl)- 
(T6). Another important result obtained in Table [ill is 
that the lowest spin triplet excitation gap decreases as the 
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FIG. 10. (Color online) Two dimer covering patterns of the super cluster for the VBS order with staggered (left panel) and 
uniform (right panel) hexagonal resonance. Super cluster (black solid lines) is connected with its surroundings through the 
empty triangle inter-dimer interaction (dashed lines) only. In the figure, thick (blue) lines denote the spin singlet pairs in the 
VBS order. 



number of topologically orthogonal inter-dimer interac- 
tions increases. Presuming that the same trend remains 
valid up to the super clusters, we expect that the spin 
triplet excitation energy gap for the super cluster gets 
smaller than that of (T6), 0.243 J, and becomes similar 
to that of the 36-site ED result with periodic boundary 
condition, 0.162 J. 

Covering entire topologically orthogonal inter-dimer 
interactions in the 36-site unit cell, the super cluster is 
connected with its surroundings through empty triangle 
inter-dimer interactions (dashed lines in Fig. 1101) only. As 
mentioned in Sec. Ill Ai the main role of the empty tri- 
angle inter-dimer interaction is to minimize the ground 
state energy. Considering that the energy-lowering effect 
of the empty triangle inter-dimer interaction has been 
rather well captured by the bond operator mean-field the- 
ory, it is reasonable to assume that the same is true for 
the spin cluster operator mean-field theory. For this rea- 
son, for future work, we propose to analyze the Kagome 
lattice antiferromagnetic Heisenberg model via the spin 
cluster operator theory with super clusters. 



VIII. SUMMARY 

In this paper, we have considered the valence bond 
solid state with a 36-site unit cell as a possible ground 
state for the antiferromagnetic Heisenberg model on the 
Kagome lattice. Our study is motivated by the real- 
ization that there is absolutely no quantum fluctuation 
about the spin-singlet product state when all inter-dimer 
interactions are topologically orthogonal. Assuming that 
the ground state tries to minimize the quantum fluctua- 
tion by forming as many topologically orthogonal inter- 
dimer interactions as possible, we arrive at the valence 
bond solid state with a 36-site unit cell. The valence bond 
solid state has been analyzed previously by the bond op- 



erator theory where quantum modes emerging from a 
bond connecting two nearest neighboring spins are re- 
garded as a natural building block for the dynamics of 
the system. While the bond operator theory is able to 
correctly reproduce the valence bond solid order in the 
ground state, it is shown that the spin triplet excitation 
spectrum is grossly overestimated in energy. 

To diagnose the problem, we have carefully investi- 
gated effects of the two inter-dimer interaction struc- 
tures, each named as the topologically orthogonal and 
the empty triangle structure. From exact diagonaliza- 
tion of various finite-sized spin clusters, it is learned that 
one has to treat the topologically orthogonal inter-dimer 
interaction in a more accurate manner than in the bond 
operator mean-field theory. As a first step toward the ac- 
curate description for the topologically orthogonal inter- 
dimer interaction, the spin cluster operator theory is de- 
veloped, in which natural quantum modes of the sys- 
tem are chosen from the spin clusters containing as many 
topologically orthogonal inter-dimer interactions as pos- 
sible while preserving the symmetry of the VBS order 
ground state. In this paper, the spin cluster operator 
theory is used to analyze two different VBS states with 
one having the staggered hexagonal resonance and the 
other the uniform counterpart. 

As a result, it is shown that the lowest spin triplet 
excitation gap is much reduced from that of the bond 
operator theory, closely approaching the numerical esti- 
mate obtained via exact diagonalization as well as the 
density matrix renormalization group method. Despite 
the better agreement, the fact that the lowest spin triplet 
excitation mostly comes from the inner region of the unit 
cell involving the pinwheel structure may be an artifact 
of the current spin cluster operator theory that uses a 
subset of the unit cell as a spin cluster. For future work, 
we propose to use super spin clusters covering the entire 
unit cell for the spin cluster operator representation 
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Appendix A: Mean field decoupling of the quartic 
term 



In this section, we describe how to decouple the quar- 
tic interaction term between spin triplet bond particles 
into a mean-field quadratic Hamiltonian. Restoring the 
summation symbol, the usual quartic term of the spin 
triplet bond operator can be written as follows: 

— ^aPf^aiiv't\i3'tl'i't\iJ'2u 

— ^ (^1,34/3^17^27 — t\pt2pti^t\^ 



^ (4^40)(*l7*27)) 



/3#7 



{t\pt2p){tlA^)) 



= ^{Q*Q-P*P) 



+ 3 E {P^Lho. - QtljL + H. c.) , (Al) 

a 

where P — {iiai\a) and Q — {tiat2a)- Note that, in 
the first step of the above equation, we have used the 
identity that ^ai3-i^a^iu = Sp^^S^y^ - Spv^^^. 



Appendix B: Bond operator theory in the empty 
triangle structure 

The Heisenberg-model Hamiltonian for the empty tri- 
angle unit structure. El, in Fig. [2]is given by 

He1,J — J (Sq • Sb + Sc • Srf + Se • Sf) 

+ XJ{Sb-Sa + Sd-Sf + SfSb). (Bl) 

Performing essentially the same procedure in the analysis 
of the topologically orthogonal structure, the bond oper- 



ator mean field Hamiltonian can be obtained as follows: 

3 



J 



1=1 

+ XJJ2 {Mj^+ic^ + Btljl+i^ + H. c.) , (B2) 



where 
cei.o = 3 



A = 



B 



-P 



(B3) 

(B4) 
(B5) 



Note that, while there are three different bonds in El, we 
set the spin-singlet condensation density to be identical 
across all the bonds with the same being true for the 
corresponding chemical potential. That is, si = S2 = 
S3 = s and fJ-i — fJ-2 = fJ-s = M- In the above, the mean 
field parameters, P and Q, are defined as follows: 



P \ tjf-yt. 



(B6) 
(B7) 

where „ is equal to tia. 

The mean field Hamiltonian, , is block- 

diagonalized in a proper basis set that takes into account 
the threefold rotational symmetry of the system. The 
new basis set of the triplet operator is defined as follows: 



(B8) 



where z„i = 3"™ and 771 = —1,0,1. In terms of this 
basis set, the mean field Hamiltonian is given by 



J^Ei — ^Ei.o — 



+ 2 E 



P t 

''ma '' — ma 



An 

Bl 



Bm 

A-r 



(B9) 



where 



Ara = ^ - M + XJ{Z-raA + Z,„A*), (BIO) 



Brn — XJ{z„i + Z^m)B. 



(Bll) 



A complete diagonalization of is performed via 

the Bogoliubov transformation: 

1 

H^f = egr,Ei + '^™7l7m , (B12) 
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where tOm is the positive eigenvalues of the following 
eigenvalue problem: 



-Bl 



,(1) 



(1) (2) 




(1) (2) 
(1) (2) 







and 



9 /J 



£gr.El — ££1,0 ~ 2 W ^ 



2 ^ 

m— — 1 



(B13) 



^m- (B14) 



Note that the Bogoliubov quasi-particle operator is trans- 
formed to the original triplet bond operator through the 
transformation matrix composed of the eigenvectors in 
the above: 




(B15) 



where Um'^'' and Um'^^ are normalized in such a way that 



'■in I 
,(2)|2 



.,(1)|2 - 



1, 



.,(2)|2 _ 



= -1, 



(B16) 
(B17) 
(B18) 

The spin-spin correlation value can be computed for 
the (a, b) and {b, d) link in the ground state as follows: 



y(l)*u(2) „ „(1)*„(2) = 



(S, ■Sb)= ( gr 

El 



3 I 1 ^ — A I 



gr 



and 



(Sb ■ ^d) — {Hei,j) 



- ^e,..Ei 



(B20) 
(B21) 



The spin-spin correlation value for the other links can be 
determined through the space rotation symmetry. 

The spin-spin correlation value for the excited state 
can be computed in a similar manner: 



El 



gr 



I ma. 



gr 



Finally, {Si, ■ Sd) can be obtained in a similar manner 
described in Eq. ((B20| : 



3J 



(egr,El 



(B24) 



Appendix C: Brief sketch for the derivation of the 
spin cluster operator representation 

The spin cluster operator representation shown in 
Eq. (|32|) is a result of the simplification due to a spin 
rotation symmetry. With the spin rotation operator de- 
fined as 

R«(<^) = e"**"■^ (CI) 
the spin triplet states satisfy the following relationships: 



. (I) 



R& (tt) \tp) = 



(C2a) 



(C2b) 



where q;,/3,7 G {a;,y, z}. Meanwhile, the spin rotation 
operator annihilates the spin singlet state, \s): 



Rrl(0)|s) 



(C3) 



Also, the spin operator itself satisfies the following trans- 
_i formation rules: 



R.ws„,Rtw = {f;; . (C4b) 

Making use of the above results, one can derive a series 
of general properties for the matrix element of the spin 
operator: 



0, 



El 



= s 

4 



^ 51 El (gr| 7™a4a'^la'7,L |gI")El ' 

a.' 

(B22) 



(C5) 

(s|^„a|i^)=0 (C6) 

(s| S'na; \tx) = (s| 'S'ny = (s| . (C7) 

For any two triplet states, \tia) and \t213), 

(^li/l Snx \t2z) — — {tlz \ Snx \t2y} 

= {tlz \ Sny \t2x) — — {tlx \ Sny \t2z) 

= {tlx \ Snz |^2i/) = ^ {tlyl Snz \t2x) ■ (C8) 

Note that 



where the last term is evaluated by using Eq. (|B8p and 
(|B15|) as well as the fact that 7™^ |gr)Er =0: 



(il^l Sna 1^27) — 



(C9) 



El (gr| 7™aiIa'^l«'7rL |gr)Ei 
1 / ... 2 

3 



(1) 



(2) 



i = -l 



,(2) 



if any two indices among (a, /3, 7) are the same. Finally, 
it is due to Eq. (jC8|) that the spin operator representa- 
tion matrices in Eq. (j35p and (|42p are anti-hermitian. 



(B23) 
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Appendix D: Details of the mean field Hamiltonian 
in Section [iTll 

The mean field Hamiltonian ([50| has following form: 

J 



Hmf 

M 

k n=l 

iV 

k n=l 

AT 

+ EE(^" 

k n=l 

M N 

+ ^-^E E E(^-('') + =«4n„(k))t^Uk)c^t(_k) + H. c. 

Ic m— 1 n— 1 

A/ 

+ ^-^E E E(^™"(k) + ^™„(k))ea(k)C'(k) +H. c. 

k m— 1 n— 1 

A/ iV 

+ ^-^E E E(C'™«(k) +^™„(k))ejk)C^(-k) +H. c. 

k m— 1 n— 1 

M N 

+ ^-^E E E(^"«(k) + %nn{k))tgi{k)t^^{k) + H. C. 
k m— 1 n— 1 

iV iV 

+ ^-^E E E(^"-(k) +,^„,„(k))t^it(k)iH2t(_k) +H. c. 

k m— 1 n— 1 

+ ^^E E E(^™«W + ^™„(k))0(k)Ci(k) + H. c, (Dl) 

k m=l n— 1 



where (e^) means the energy of triplet state \tna) in the David star (extended hexagon), and Eo is defined in (|46p . 
Remaining undefined quantities are defined below: 
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M N 



k,l=l p,q=l 
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r r 

P 9 



M N 
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M N 
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4-4 -l\ 
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-1 



ZkZi 



2 — 2I 
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M AT 
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M N 
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M N 
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M JV 



^™.(k) = -^EE4'«'+^^^^) 
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(D2) 

(D3a) 
(D3b) 
(D3c) 
(D3d) 
(D3e) 

(D3f) 



(D4a) 



(D4b) 



(D4c) 



^,^^2(0)] , (D4d) 
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k=l 1=1 

= E E + ^l^'-n ^) + + &!^[-™-.T^)" 

A:=l (=1 

• [rfr-'.^^,"{-rc)e^^-^- + r^r,-i^,5«(r^)e-'^"-- + rlr^'nr-' ^l^^" {0)] . (D4f) 
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